对流扩散方程的QC3无网格法

对流扩散方程作为偏微分运动方程的分支,在流体力学、气体动力学等领域有着重要应用.为解决对流扩散方程难以通过解析法得到解析解的难题,采用二阶一致3点积分(Quadratically Consistent 3-Point Integration,简称QC3)提高无网格法的计算效率,通过对积分点上形函数导数的修正,改善无网格法的精度和收敛性.本文将QC3无网格法拓展到对流扩散方程问题中,时域离散采用广义特征线Galerkin法,空间离散采用QC3法.数值结果表明,应用QC3无网格法得到的对流扩散问题数值解与解析解十分接近,验证了QC3无网格法解决对流扩散问题的可行性.     

移动加热器法生长碲锌镉晶体的组分输运与界面形貌研究

移动加热器法(THM)生长碲锌镉晶体时,界面稳定性对晶体生长的质量有很大影响。本文基于多物理场有限元仿真软件Comsol建立了THM生长碲锌镉晶体的数值模拟模型,讨论了Te边界层与组分过冷区之间的关系,对不同生长阶段的物理场、Te边界层与组分过冷区进行仿真研究,最后讨论了微重力对物理场分布的影响,并对比了微重力与正常重力下的生长界面形貌。模拟结果表明,Te边界层与组分过冷区的分布趋势是一致的,在不同生长阶段,流场中次生涡旋的位置会发生移动,从而导致生长界面的形貌随着生长的进行发生变化,同时微重力条件下形成的生长界面形貌最有利于单晶生长。因此,在晶体生长的中前期,对次生涡旋位置的控制和对组分过冷的削弱,是THM生长高质量晶体的有效方案。     

DFT modeling of indoline spiropyrans with a cationic substituent in the gas phase

Isomeric forms of indoline spiropyrans show unusual behavior compared with similar compounds, according to experimental data. DFT modeling for gas phase was made to consider the simplest case without environmental effects, which revealed the intramolecular reasons for occurrence of ring opening reaction depending on the particular structure of the compound. The questions of charge redistributions, the changes of geometry and chemical bonds in the structures are also discussed.     

The Decisive Role of Spin States and Spin Coupling in Dictating Selective O2 Adsorption in Chromium(II) Metal–Organic Frameworks**

The coordinatively unsaturated chromium(II)-based Cr₃[(Cr₄Cl)₃(BTT)₈]₂ (Cr−BTT; BTT³⁻=1,3,5-benzenetristetrazolate) metal–organic framework (MOF) has been shown to exhibit exceptional selectivity towards adsorption of O₂ over N₂/H₂. Using periodic density functional theory (DFT) calculations, we attempted to decipher the origin of this puzzling selectivity. By computing and analyzing the magnetic exchange coupling, binding energies, the partial density of states (pDOS), and adsorption isotherms for the pristine and gas-bound MOFs [(Cr₄(X)₄Cl)₃(BTT)₈]³⁻ (X=O₂, N₂, and H₂), we unequivocally established the role of spin states and spin coupling in controlling the gas selectivity. The computed geometries and gas adsorption isotherms are consistent with the earlier experiments. The binding of O₂ to the MOF follows an electron-transfer mechanism resulting in a Cr^III superoxo species (O₂^.−) with a very strong antiferromagnetic coupling between the two centers, whereas N₂/H₂ are found to weakly interact with the metal center and hence only slightly perturb the associated coupling constants. Although the gas-bound and unbound MOFs have an S=0 ground state (GS), the nature of spin the configurations and the associated magnetic exchanges are dramatically different. The binding energy and the number of oxygen molecules that can favorably bind to the Cr center were found to vary with respect to the spin state, with a significant energy margin (47.6 kJ mol⁻¹). This study offers a hitherto unknown strategy of using spin state/spin couplings to control gas adsorption selectivity in MOFs.     

Effect of wall slip on the viscoelastic particle ordering in a microfluidic channel

The formation of a line of equally spaced particles at the centerline of a microchannel, referred as “particle ordering,” is desired in several microfluidic applications. Recent experiments and simulations highlighted the capability of viscoelastic fluids to form a row of particles characterized by a preferential spacing. When dealing with non-Newtonian fluids in microfluidics, the adherence condition of the liquid at the channel wall may be violated and the liquid can slip over the surface, possibly affecting the ordering efficiency. In this work, we investigate the effect of wall slip on the ordering of particles suspended in a viscoelastic liquid by numerical simulations. The dynamics of a triplet of particles in an infinite cylindrical channel is first addressed by solving the fluid and particle governing equations. The relative velocities computed for the three-particle system are used to predict the dynamics of a train of particles flowing in a long microchannel. The distributions of the interparticle spacing evaluated at different slip coefficients, linear particle concentrations, and distances from the channel inlet show that wall slip slows down the self-assembly mechanism. For strong slipping surfaces, no significant change of the initial microstructure is observed at low particle concentrations, whereas strings of particles in contact form at higher concentrations. The detrimental effect of wall slip on viscoelastic ordering suggests care when designing microdevices, especially in case of hydrophobic surfaces that may enhance the slipping phenomenon.     

The Role of Radical Bridges in Polynuclear Single-Molecule Magnets

Employing radical bridges between anisotropic metal ions has been a viable route to achieve high-performance single-molecule magnets (SMMs). While the bridges have been mainly considered for their ability to promote exchange interactions, the crystal-field effect arising from them has not been taken into account explicitly. This lack of consideration may distort the understanding and limit the development of the entire family. To shed light on this aspect, herein we report a theoretical investigation of a series of N -radical-bridged diterbium complexes. It is found that while promoting strong exchange coupling between the terbium ions, the N -radical induces a crystal field that interferes destructively with that of the outer ligands, and thus reduces the overall SMM behavior. Based on the theoretical results, we conclude that the SMM behavior in this series could be further maximized if the crystal field of the outer ligands is designed to be collinear with that of the radical bridge. This conclusion can be generalized to all exchange-coupled SMMs.     

Catalytic Nitrene Transfer by an FeIV-Imido Complex Generated by a Comproportionation Process

Nitrene transfer reactions have emerged as one of the most powerful and versatile ways to insert an amine function to various kinds of hydrocarbon substrates. However, the mechanisms of nitrene generation have not been studied in depth albeit their formation is taken for granted in most cases without definitive evidence of their occurrence. In the present work, we compare the generation of tosylimido iron species and NTs transfer from Fe^II and Fe^III precursors where the metal is embedded in a tetracarbene macrocycle. Catalytic nitrene transfer to reference substrates (thioanisole, styrene, ethylbenzene and cyclohexane) revealed that the same active species was at play, irrespective of the ferrous versus ferric nature of the precursor. Through combination of spectroscopic (UV-visible, Mössbauer), ESI-MS and DFT studies, an Fe^IV tosylimido species was identified as the catalytically active species and was characterized spectroscopically and computationally. Whereas its formation from the Fe^II precursor was expected by a two-electron oxidative addition, its formation from an Fe^III precursor was unprecedented. Thanks to a combination of spectroscopic (UV-visible, EPR, Hyscore and Mössbauer), ESI-MS and DFT studies, we found that, when starting from the Fe^III precursor, an Fe^III tosyliodinane adduct was formed and decomposed into an Fe^V tosylimido species which generated the catalytically active Fe^IV tosylimide through a comproportionation process with the Fe^III precursor.     

Mechanism of BPh3-Catalyzed N-Methylation of Amines with CO2 and Phenylsilane: Cooperative Activation of Hydrosilane

BPh₃ catalyzes the N-methylation of secondary amines and the C-methylenation (methylene-bridge formation between aromatic rings) of N,N-dimethylanilines or 1-methylindoles in the presence of CO₂ and PhSiH₃; these reactions proceed at 30–40 °C under solvent-free conditions. In contrast, B(C₆F₅)₃ shows little or no activity. ¹¹B NMR spectra suggested the generation of [HBPh₃]⁻. The detailed mechanism of the BPh₃-catalyzed N-methylation of N-methylaniline ( 1 ) with CO₂ and PhSiH₃ was studied by using DFT calculations. BPh₃ promotes the conversion of two substrates (N-methylaniline and CO₂) into a zwitterionic carbamate to give three-component species [Ph(Me)(H)N⁺CO₂⁻⋅⋅⋅BPh₃]. The carbamate and BPh₃ act as the nucleophile and Lewis acid, respectively, for the activation of PhSiH₃ to generate [HBPh₃]⁻, which is used to produce key CO₂-derived species, such as silyl formate and bis(silyl)acetal, essential for the N-methylation of 1 . DFT calculations also suggested other mechanisms involving water for the generation of [HBPh₃]⁻ species.     

Shedding vortex simulation method based on viscous compensation technology research

Owing to the influence of the viscosity of the flow field, the strength of the shedding vortex decreases gradually in the process of backward propagation. Large-scale vortexes constantly break up, forming smaller vortexes. In engineering, when numerical simulation of vortex evolution process is carried out, a large grid is needed to be arranged in the area of outflow field far from the boundary layer in order to ensure the calculation efficiency. As a result, small scale vortexes at the far end of the flow field cannot be captured by the sparse grid in this region, resulting in the dissipation or even disappearance of vortexes. In this paper, the effect of grid scale is quantified and compared with the viscous effect through theoretical derivation. The theoretical relationship between the mesh viscosity and the original viscosity of the flow field is established, and the viscosity term in the turbulence model is modified. This method proves to be able to effectively improve the intensity of small-scale shedding vortexes at the far end of the flow field under the condition of sparse grid. The error between the simulation results and the results obtained by using fine mesh is greatly reduced, the calculation time is shortened, and the high-precision and efficient simulation of the flow field is realized.     

改进型单晶炉提拉系统偏心平衡研究

随着光伏行业的快速发展, 对硅单晶的品质和长晶装备的稳定性的要求也不断提高。直拉法是生产硅单晶的主要方法,通过提高单晶炉副室的高度以扩大单晶硅的生产规模。由于副室高度的大幅增加,且单晶炉提拉头质心相对于旋转轴心有一定距离,对单晶炉整体稳定性有较大影响,从而降低了单晶硅的生产质量。针对此问题,对单晶炉建立可靠的力学分析模型,采用数值仿真方法,对单晶炉整体进行动力学响应分析,计算得到副室高度增加后的单晶炉工作时中钨丝绳下端晶棒的运动规律以及最大摆动幅度,为改进设计提供依据。数值仿真分析表明提高单晶炉副室高度后,提拉头较大的质心偏心是单晶炉提拉系统发生摆动的主要原因。在此基础上提出在提拉头上添加质心调节装置,通过控制系统调节可保证提拉头质心位置在旋转轴线上以降低提拉系统的摆动。